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ABSTRACT 

We study coupled axial and polar axisymmetric oscillations of a neutron star endowed 
with a strong magnetic field, having both poloidal and toroidal components. The 
toroidal component of the magnetic field is driving the couphng between the polar and 
axial oscihations. The star is composed of a fluid core as weh as a sohd crust. Using 
a two dimensional general relativistic simulation and a magnetic field B — lO^^G, we 
study the change in the polar and axial spectrum caused by the couphng. We find 
that the axial spectrum sujfers a dramatic change in its nature, losing its continuum 
character. In fact, we find that only the 'edges' of the continua survive, generating 
a discrete spectrum. As a consequence the crustal frequencies, that in our previous 
simulation could be absorbed by the continua, if they were embedded inside it, are 
now long living oscillations. They may lose their energy only in the very special case 
that they are in resonance with the 'edges' of the continua. 
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1 INTRODUCTION 



In the last few years, a great attention was devoted to the interpretation of the quasi periodic oscillations (QPOs) discovered 
in the tail of giant flares in Soft Gamma Repeaters (SGRs). Those oscillations cover a wide range of frequencies going from a 
few Hz up to kHz. Until now only in two SGRs it was possible to observe clearly QPOs: in the SGR 1806-20 and in the SGR 
1900+14 (see |Israel et al.| ( |2005| , |Watts fc Strohmayer| ( |2006[ ) for a review). In the SGR 1900+14 four frequencies have been 
detected: 28, 53, 84, and 155 Hz while in the SGR 1806-20 several QPOs have been discovered: 18, 26, 30, 92, 150, 625 and 



1840 Hz. A recent reanalysis of the data for the SGR 1806-20 carried on by Hambaryan et al. (2011 1 with a different method. 



shows the presence of additional QPOs at 16.9, 21.4, 36.4, 59.0 and 116.3 Hz. 

Understanding the nature of the QPOs will offer an unique opportunity to investigate the properties of neutron stars 
and to constrain the nuclear matter in neutron stars as well as the topology and strength of their magnetic field. For those 
reasons, since their discoveries, QPOs have been the subject of an intensive study. The first hypothesis, in order to explain 
QPOs, was to consider pure shear modes of the crust, excited after the burst generated by giant flares. This idea proposed 
by Duncan (19981 was investigate in several analytical and numerical works (see for example Sotani et al. (20071, Samuelsson 



& Andersson (2007l). Those works evidenced that not all the lower QPOs could be explained as pure crustal frequencies. 



An alternative scenario has been then proposed by Levin (20061 and Glampedakis et al. (20061, involving giobal crust-core 
oscillations of the star. In particular, |Levin| ( |2007[ ) showed with a toy model that the QPOs spectrum should be continuum, 
with turning points that were called 'edges' of the continuum. FoUowing this idea, a significant amount of work has been 



devoted in understanding global Alfven oscillations in neutron star, using general relativistic models (Sotani et al. (20081 



Colaiuda et al. (20091, Cerdà-Duràn et al. (2009l) or Newtonian model (Lee (2008l). All those works made clear that global 
Alfven modes in a pure fiuid neutron star could not explain the small gap between the lower frequencies observed in QPOs. 

A further step was then needed: to add a solid crust and to study the effect of its presence on the oscillations of the 
core ( see |Kokkotas et al.| ( |2010[ )). |van Hoven fc Levin| ( |2011[ ), using a non-relativistic model, showed that the presence of 
the crust strongly influences the continuum spectrum. In particular, they found that near the edges of the continua other 
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discrete Alfven modes appear. This results was partially confirmed by Gabler et al. (2011 1 and Gabler et al. (2012 I, using a 
non-linear relativistic method: the authors could find some frequencies in the gaps between two near continua, however they 
stated that those oscillations were quickly damped. In a hnear relativistic simulation, Colaiuda & Kokkotas (2011 ) confirmed 
the presence of those discrete Alfven modes in the gap between two continua and showed that the crustal modes could hve 
long if they are located in those gaps. The discrete Alfven modes as well as the crustal modes in the gap between the first two 
continua could explain the somehow dense spectrum of the observed lower QPOs. However as pointed out in |van Hoven fc| 



Levin (20121, it was still difficult to explain the higher frequencies as the 625Hz because at higher frequencies the continuum 



dominates, absorbing all the discrete frequencies, and in addition the continua overlap to each other: in this way it is not 
possible to isolate the edges of the continuum and to have a unique interpretation for the 625Hz QPO. 

A possibihty to explain the high frequency QPOs came from the study of polar oscillations in magnetars. The work by 
|Sotani fc Kokkotas| ( |2009l ) showed that the polar Alfven spectrum has a discrete nature and, in addition, the lower polar Alfven 
modes have frequencies around few hundred Hz. However, it was not clear how this result could affect the axial spectrum. 

In this work, we make an attempt to investigate the effect of the coupling between axial (torsional) and polar oscillations 
on the spectrum. Our magnetar model is composed of a fluid core and a solid crust and is permeated by a magnetic field with 
strength B = 10^*^0 which acquires both poloidal and toroidal components. The toroidal component is the coupling 'pipeline' 
between the axial and polar part of the spectrum. By using a two dimensional relativistic, linear code, we find that the axial 
spectrum is strongly modified by the presence of the coupling: its continuum feature is destroyed, leaving behind only its 
'edges'. In this new picture, both Alfven and crustal modes have a discrete nature and the latter ones, if out of resonance 
with the 'edges' of the continuum, could live long enough to be seen. 

The paper is organized as follows: in Section 2 we present our background model, while in Section 3 we derive the 
perturbation equations and test our code by comparing it with earlier results. In Section 4 we study a configuration where the 
coupling between axial and polar oscillations is active only on the axial part (i.e. we do not take into account the coupling to 
the polar part). In Section 5 the fuU coupling between axial and polar oscillations is presented, including boundary conditions 
and initial data. The Section 6 and the Section 7 are devoted to comments and conclusions. 



2 BACKGROUND MODEL 

We consider a non-rotating spherical symmetric star, whose metric is given by: 



ds 



-2* 7,2 . 2A 7 2 . 2/ 7fl2 . . 2 

-e dt + e dr + r {dti + sm 



'■) 



(1) 



We restrict our model to axisymmetric perturbations: all physical quantities are not dependent on the (j> coordinate. The only 
non-zero term of the unperturbed four-velocity is m*, i.e. it^ = (e-*, 0,0,0). Finally, we do not consider the deformations 
induced on the star by the presence of a strong magnetic field, because the magnetic field has an energy Em which is a few 
orders of magnitude smaller than the gravitational binding energy Eg. Typically, Em/Eg ~ 10~''(-B/10^''G)^ (see Colaiuda 



et al. ( 2008 1 and Haskell et al. ( 2008 1 for numerical results about magnetars deformations) . 

We assume that the star is constituted of a perfect fluid in the ideal MHD approximation and of a thin crust, that, as 
flrst approximation, can be described by a shear velocity vs: 

(2) 

where e is the density and fi is the shear modulus. Note that the presence of nonuniform nuclear structure (the so called 



vs = — = i-U cm s 
e 



'pasta-phase') in the crust could reduce the frequencies of the shear modes, as it was recently suggested by Sotani (2011 1 



In this work, we neglect the superfluid component present both in the crust and in the core as well as the possible presence 



of proton superconductivity in the outer core (see Glampedakis et al. (2011 1), although we know that their appereance could 



strongly influence the features of the oscillations spectrum both in a qualitative way, changing its nature from continuum to 



discrete, and in a quantitative way, shifting the frequencies (for more details see Andersson et al. (2009 1). 
The magnetic fleld inside the star is described by the Grad-Shafranov equation: 

2 



e -f ($ — A j e 



ai 



^2 -2* 



ai = —4nji 



(3) 



where ji and ai are respectively the radial component of the four-current and the radial component of the magnetic fleld. 
Solving the equation ([3| with the appropriate boundary conditions, we get the magnetic fleld conflguration inside the star. 
In our model, the magnetic field has both poloidal and toroidal components and it is given by: 



H'' = 



/Tvr^ 



-ai, 



sin eC 



-ffli. 



(4) 
(5) 



where is a parameter that describes the ratio of the toroidal magnetic fleld respect to the poloidal one, see |CoIaiuda et aLj 
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Figure 1. The eigenfunction of the fundamental polar mode for the function w(t,r,6) at 300Hz for £ = 2 on the left and of the 
fundamcntal polar mode at 500Hz for £ = 4 on the right. Those plots arc in good agrccmcnt, rcspcctivcly, with thc first plot and thc 
third plot in figure 3 in l |Sotani &: Kokkotas|2009l l. 



(20081. The interior magnetic field is matched with a pure poloidal magnetic field outside the star. This choice implies that 



the toroidal magnetic field must be set equal to zero outside the star: in this way, a surface current is created. 



3 PERTURBATION EQUATIONS 

We consider both polar and axial perturbations, using the Cowling approximation, i.e Sg^^ = so the components of the 
four-velocity u'^ are given by |Sotani fc Kokkotas| pOOOt : 



Su^ 
ou = — r 
Su''' = e 



0, 



r-^e-'^-'^dtW{t,r,e), 



-^e-'^dtV{t,r,e), 
■*-''9tJ(t,r,e). 



(6) 
(7) 
(8) 
(9) 



The polar perturbations involve perturbations of the density and the pressure. For adiabatic perturbations and using the 
first law of the thermodynamics, the energy density and pressure variation could be written as: 



ot = (e + p) dre^e 

n r-' 



&p - 



An 



IP- 



OrP^e 



(10) 
(11) 



where 7 is the adiabatic constant and An/n is the Lagrangian perturbations of the baryon number density. 



With these assumptions, the Hnearized equations of motion are given by Sotani et al. ( 2008 1 
(e + p + H'^)5u''.,u'' = {5e + 5p + ^SH^Hc) u'^.Ju" 



+ {u'^SUa + S'^Ua) 



H'^H" - 



p^\h^ 



+ K 



H'^SH" + 5H'^H" - g''"{5p + Hp5H^) 



(12) 



where /i^j^ = g^j^ + UpU^ and ST'^'f^ is the Hnearized shear stress tensor (see the Appendix for the fuU derivation of the 



shear stress tensor). The equation (12l leads to three partial differential equations (PDEs) for /i = r,e,(j), which have to be 
solved with the appropriate boundary and initial conditions. 



3.1 Comparison with earlier results 

In order to check our code for the polar oscillations, we implement the equation (22) of [Sotani fc Kokkotas| ( [2b09[ ) considering 
an axisymmetric non-rotating neutron star, permeated with a pure poloidal magnetic field. To be consistent with the model 
used in Sotani fc Kokkotas ( 2009 1, the star is constituted of a pure fluid without the presence of a soHd crust. To get the code 
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stable, we use an artificial viscosity of the second order. Note that in Sotani & Kokkotas (20091 the artificial viscosity is of 



the fourth order. Despite this difference, our results are in good agreement with the ones found by Sotani & Kokkotas ( 2009 1 
as it can be seen from figure [l] where the eigenfunctions for the polar modes for 1 = 2 and £ = 4 are plotted. 



4 FIRST CONFIGURATION 



In order to investigate the couphng between polar and axial osciUations: we should solve the system of equations described 
in | |12[ ) which involve both the polar [fi = r, 9) and axial part (/! = </!>), obtaining three coupled but different PDEs for the 
functions 3^, W and V . Each of those equations has the form: 



attF(i, r, 0) = V F + coupling terms (13) 

where F could be one of the function y, W and V and the coupHng terms involve the derivative both in space and time of the 
remaining two functions. The equations that we get are quite lengthy and require a significant amount of computational work 
to be solved. In order to speed up our computational time and as first attempt to understand the effect of the coupHng, we 



choose to evolve the full equation ( 13 1 only for the axial function y, while for V and W we solve the equations ( 13 1 without 



the presence of the coupling terms, setting them artificiaUy equal to zero. 



In this approximation, the equation (12| reduces for /i = r,0 to the one found in Sotani & Kokkotas (20091. However, 



in addition to Sotani & Kokkotas ( 2009 1 we consider the presence of a sohd crust while the magnetic field has not only a 



poloidal component but also a toroidal one. As we explain later, this is needed in order to provide a couphng between the 
polar and axial osciUations: notice that without the presence of toroidal field, there is no couphng between the classes of 
perturbations. An aUernative way to have a consistent couphng between the polar and the axial oscillations is to consider 
the deformation that the magnetic field induces on the shape of the star: however, this means to consider a perturbation 
of the metric given in ([l]), that wiU lead to a much more comphcated structure of the final equation to solve and probably 
considerable weaker coupUng. Furthermore, there is a general belief that a combined poloidal-toroidal magnetic field could 



form a stable configuration (see Braithwaite & Nordlund (2006l) 
The analytic form of equation ( 12 1 for fi = (j) is: 

d^y 
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In this equation the coefficients Aoo, A20 
given by 



2-Kr'^ 
Aii, A02 



2t7 
ai V,r 



Aoi, A\o depend on the coordinates r and but not on time. They are 
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(18) 
(19) 
(20) 



As we aheady pointed out, the polar contribution in equation ( 14 1 is always coupled with the toroidal field: in the absence 



of a toroidal field then the equation reduce to the one used for the study of pure torsional osciUations ( [Sotani et al.|[2007l ). 
Then equation ( 14 1 must be solved together with equation (22) in Sotani & Kokkotas ( 2009 1 for the polar functions W and 



V. However, in Sotani & Kokkotas (20091 the star consists of a pure fluid without crust: this we have to calculate the crust 
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contribution for the polar perturbations (see the Appendix for the fuU derivation) . The final forrn of the polar perturbation 
equations to be solved is: 



- {5e + 5p) $'e-^^ + H''5H"^ 
+ [-2<i>'e-^''Ho + H\o 
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(22) 



The equations above must be solved with the appropriate boundary conditions hsted in section [4~T| 

We want to point out that this type of configuration aUows the transport of energy from the polar osciUations to the 
axial osciUations but not vice-versa. We wiU then expect to find the trace of the polar spectrum in the axial osciUations but 
no actual change in the polar spectrum. This is indeed what we find: the axial spectrum presents a significant excitation in 



the frequencies correspondent to the polar ones found in Sotani & Kokkotas (20091. Those frequencies do not present any 



significant excitation in the case of pure torsional osciUations, see Colaiuda & Kokkotas (20111. On the contrary, the polar 
spectrum does not show any change, as expected. 

This study demonstrates that, in order to understand how the polar oscillations are affected by the coupUng to the axial 



ones, we need to study the fuU problem, i.e. consider the equation ( 13 1 also for the polar functions W and V 



4.1 Boundary conditions 



The equation ( 14 1 must be solved with the appropriate boundary conditions for the axial function y as well as for the polar 



functions 'W and V. In order to simpUfy the boundary conditions on the surface, we use two new function w(t,r,e) and 



v{t,r,e) defined as (see Sotani & Kokkotas (20091 



w = eW V — eV 
With this substitution, the equation becomes: 
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Using those new functions the boundary conditions are chosen as foUowing: 
at the centre we require the regularity of the equations: 

:y^(r = 0)=0, ii;(r = 0) = 0, ^(r = 0) = , 
at the surface, we use the zero traction condition: 

dry = 0, drv{r = i?) = 0, drw{r = i?) = , 
on the magnetic axis at = 0, we require axisymmetry: 

dey = 0, dew = 0, dev = . 



v,e 



(24) 



(25) 
(26) 
(27) 
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• At the equatorial plane, we require equatorial plane symmetry a,t 9 — 7r/2: 

dgy = 0, dgw = 0, dgv = 0, 

or antisymmetry: 

y = o, w = 0, v = 0, 

• at the crust-core interface, we require the continuity of the functions: 



dry = 



1 + 



VA 



dry. 



(28) 



(29) 



(30) 



drV 



1 + 



va 



drV , 



(31) 



drW = 



1 + 



VA 



drW . 



(32) 



where Vs, the shear velocity, and va, the Alfven velocity, are respectively given by: 

Vs = [ -] ~ 10^^ cm/s, VA = ^ 



(33) 



where n is the shear modulus. 



5 POLAR AND AXIAL OSCILLATIONS: COUPLING 



We now derive explicitly the equations for the polar function W and V, including the presence of the crust presented in 
equations (62l and (631 and the coupling with the axial function y. Substituting the equations ([7| and ([8| in equations (21l 
and ( 22 I we get: 



dttw{AooAii ~ AioAoi) = FrAn - F^Àoi + Ay,^^^dtty + ^hoi^s^ + -4^109^:^^ 



(34) 



dttv{AooAxi - AioAoi) = F^Aoo - F^Aw + Ay^^^dtty + Ay^^dgy + Ay^^dry 
where all the coefScients are function only of the coordinates r, 9: 



(35) 



A~oo= {t + p + H''Hg + Ht'H^)r-\-^'^-'' , (36) 

Àoi = r-^H^'He , (37) 

Aio = -r-^H^Hr , (38) 

Àii= -{e + p + H''Hr + H't'H^)r-^e-^'^ , (39) 

= -Caicoter-^e-'^^Aiia'i - A)iai), (40) 

^aoo = -(^axcot9r-'^e-'^{Aooai - Axoai) , (41) 

A'y„^= 2Ca?sinċ'cos6le~'^Aii , (42) 

^aoi = 4Ca?sin6'cos6lr-2e-^'^yloo , (43) 

Ayio = 2Ca?sin6'cos6le-'^ALi , (44) 

Ay^g = 4:Calsin9cos9r-'^e-^'^Aoo, (45) 
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while the two expressions J-V and J-V are (Sotani & Kokkotas (2009l): 

Tw = - {5e + 5p) ^'e-^^ - e~^^Sp^r + f A' + H^'SH'' + (5H\, + SH^g + SH"^ .) + H'^SH^g + /f^^ff" 
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(47) 



Note that, as for the previous configuration, the coupling between the axial and the polar functions is always determined 
by the toroidal component of the magnetic field. 



5.1 Numerical method 



We adopted the same equation of state (EOS) used in Sotani & Kokkotas (20091, i.e a polytropic EOS that can reproduce 
quite well the tabulated data for the realistic EOS A ( [Pandharipande 19711. For more details about this EOS see section 
4.1 in Sotani & Kokkotas (20091. We focus on a stellar model with a mass M — 1.4AfQ and a radius R = 10.35 km while 
the value for the magnetic field strength that we assume is = 10^'' Gauss on the pole. This magnetic field is stronger 
than the magnetic field that we used in our previous simulations. However, as was already pointed out by |Sotani fc Kokkotas| 
(20091, a magnetic field of this strength is needed in order to get a cleaner imprint of the polar frequencies. We will explain 



the implications of such strong magnetic field on the spectrum later. 

In all the simulations that we present the parameter ( is set to ^ = 0.1. We perform also a simulation witli ( — 0.2 and 
we find that the energy transfer between the polar and the axial system becomes, as expected, more efficient. 

We used a numerical grid in the plane {r,e) 150 x 50, varying e from to n/2. The accuracy of the code was tested, 
performing a simulation with a 200 x 50 grid: the results show that the frequencies are not influenced significantly by a change 
in the number of grid points. The stability of the scheme was also checked: our simulation lasted few seconds and the use of 
artificial viscosity does not damp significantly the oscillations. The base of crust ( rcrust ) is taken at e = 2.4 x lO'^^ g/cm? , 



according to the crust model by Negele & Vautherin ( 1973 1 (in the foUowing NV). Core and crust are coupled via the interface 



conditions (30l-(32l 



As initial conditions, considering the requirement of regularity at the centre of the star, we use: 
wir^e) = e/+^sin6iP£(cos6i), 



v{r,e) = 

y{r,e) = 



— -r smedePe{cose) . 
er^+^ sin ePi{cos9) . 



(48) 
(49) 

(50) 



6 RESULTS & COMMENTS 
6.1 On the nature of the spectra 

We are interested to understand not only how the perturbations propagate in the star but also to determine the changes that 
the coupling between polar and axial oscillations brings on the stellar spectrum. In particular, we want to investigate if the 
nature of the spectrum (discrete for polar modes and continuous for axial ones) is preserved or if it changes substantially. In 
addition, we want to study how efficiently the oscillations are propagating when just one type of perturbation is excited. 

For this reason, we performed different runs, exciting initially just the polar oscillations and looking for the propagation 
of the perturbation on the axial part. We notice that the axial part begins to oscillate immediately after the onset of the 
perturbation. This could be easily understood if we consider what is the frequencies range of both polar and axial oscillations: 
the polar oscillations have frequencies going from 300 Hz to kHz, while the lower axial frequency has been found around 
16 — 20 Hz, depending on the EOS used. This means that, while the axial oscillations need a longer evolution in order to be 
clearly detected, the polar oscillations could appear after a much smaller amount of evolution time. However, the coupling 
between the two systems is so efficient, even for relatively small values of the toroidal magnetic field, that the energy can 
fiow quickly from one system to the other. In this way, the axial oscillations have been excited significantly even if the initial 
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50 100 150 200 250 300 350 400 

/ (Hz) 

Figure 2. The FFT of the function y taken at a fixed radius and varying the angular position of the observing point, i.e. the valuc of 
6. It is cvident that the peaks are independent from the observer's position. In addition, we can see a particular excitation in the axial 
spcctrum at high frequencies (around 200Hz), probably due to the coupling with the polar spectrum. 



perturbation was only on the polar part of the equations. However, the axial oscillations have been influenced by the presence 
of the couphng to the polar ones in a very dramatic way. 



We know from previous works (Levin (20061, Glampedakis et al. (20061, Sotani et al. (20081, Colaiuda et al. (20091 



|Cerdà-Duràn et al.|p009| , |van Hoven fc Levin| ( |2011[ ), |Gabler et al.| ( |2011| ), |Colaiuda fc Kokkotas| ( |2011| ), |Gabler eraLl ( |2012[ ) 

) that the spectrum of the axial oscillations has a continuum component generated by the presence of a strong magnetic 
field inside the star. The presence of the couphng between axial and polar oscillations changes this picture and the continuum 
reduces to a discrete spectrum. We observe this change in the foUowing way: we choose a point with certain (r, 9) coordinates 
inside the star, then we take the FFT of this point, getting a certain frequency that corresponds to a certain mode. Then, 
by fixing the value of r, we change the value of 6, i.e. we choose a different phase, and we perform again a FFT of the point, 
getting a new frequency, and so on, scanning all the 9 range. As it can be seen in figure[2] the frequencies that we get fixing 
r and varying the value of 6 do not change, confirming their discrete nature. However, comparing those new results with the 
ones that we found in Colaiuda fc Kokkotas (2011 1, we noticed that those 'new' discrete frequencies correspond to the edges 
of the continuum that we found previously. 

It is then important to examine if the new nature of the torsional oscillations spectrum is consequence of the interface 
condition between crust and core and of the coupling of the polar and axial oscillation or it is a consequence of the coupling 
alone. In order to understand this, we solve the equations (14l, (35l and (34l for a pure fluid star without a solid crust. We 
find that the new discrete character of the torsional oscillations is preserved even in this alternative scenario, showing that 
the discrete nature of the spectrum is a consequence of the coupling between polar and axial oscillations. Since the coupling 
happens through the toroidal component of the magnetic field, it can be stated that the magnetic field configuratton has a 
remarkable role in determining the final nature of the spectrum of the torsional osctllation. In fact, we find that even for small 
values of the toroidal magnetic field, the coupling between axial and polar oscillations is already effective enough to change 
the discrete nature of the torsional oscillations. 

In addition, the axial spectrum presents strong excitations at high frequencies, probably due to the coupling with the 
polar spectrum: in fact, as we already pointed out, the first lower polar frequency is located at few hundred Herz. The strong 
excitation at high frequencies for the function y can be seen in figure|2] This excitation was not present in the case of pure 
axial modes (Colaiuda fc Kokkotas (2011)). 
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Figure 3. The FFT of the polar function w{t, r, 9) (red line) and the corresponding density perturbations 5e (green line). The amplitude 
of thc density perturbation is several order of magnitude smaller than the one of thc polar function w{t, r9): this implies that the density 
perturbations are not excited to a significant level and, as consequence, the possibility of detecting gravitational waves is quite low. 



6.2 Gravitational waves 



The polar spectrum present variations too. In particular, we can see the appearance of lower frequencies that are imprints of 
the coupling with the axial osciUations. These frequencies correspond to the discrete ones of the axial spectrum. One could 
think that this last feature of the polar spectrum could have an important effect in the possibihty of detecting gravitational 



waves from magnetars. In fact, we know from the works by Levin & van Hoven (20111, Lasky et al. (20111, Ciolfi et al 



(20111, Zink et al. (20121 that the possibility of detecting gravitational waves from torsional oscillations or hydromagnetic 



instability in magnetars is quite low. However, if density perturbations are excited in the polar spectrum by the coupHng 
with the torsional one, then the possibility of observing gravitational waves could increase significantly. For this reason, we 
study the density perturbations that have been generated by the coupUng between the axial and polar oscillations through 
the fromula: 

(p + e) / W,r e,rW 1 Vfi cot9v\ 1 



(51) 



The results are shown in figurejS] the density perturbations are several orders of magnitude lower than the normal axial 
perturbations. We then conclude that density perturbations are not excited in a significant level and the possibility of detecting 
gravitational waves remains quite low. 



6.3 On the nature of some special modes 



It is worth to notice that in Colaiuda & Kokkotas (20111, two families of discrete frequencies have been found: we called 



them 'crustal modes' and 'discrete Alfven modes': while the former modes are shear axial modes that can be recovered even 
in absence of magnetic field, the latter are a consequence of the presence of the solid crust and of its coupling with a fluid 



core and they are directly proportional to the magnetic field strength. In addition, in Colaiuda & Kokkotas (2011 1, we found 



that if a crustal mode is embedded in the continuum, then its energy will be quickly damped and eventually absorbed by the 
continuum. In the configuration that we study in this paper the continuum disappears and then the crustal modes can be 
effectively damped only if they are in resonance with the discrete frequencies of the continuum. The speed of the damping 
process depends on the strength of the magnetic field: for high values of the magnetic field, the crustal frequencies, if in 
resonance with the discrete modes of the continuum, are absorbed considerably faster. 

Regarding the 'discrete Alfven frequencies', it is difficult, in this new model, to distinguish between the edges of the 
continuum and the 'discrete Alfven frequecies' since both of them have now a discrete nature and both of them are proportional 
to the magnetic field. 

Concerning the identification of the QPOs, we try to solve the puzzle of the 625 Hz QPO that has been reported in 



Strohmayer & Watts (20051, Watts & Strohmayer (20071. In our previous work we could not explain this long living QPO 



because, at high frequencies, the edges of the continua present in absence of coupling could not be isolated anymore since the 
various continua start to overlap with each other, and the crustal mode are quickly absorbed from the continua. In this new 
model, the torsional spectrum is discrete but the spectrum of the oscillations at high frequencies are too dense to allow an 
univocal explanation for the 625 Hz QPO. This QPO could still be explained as a higher polar mode. 
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r CONCLUSION 



In this paper we study the couphng between the axial and the polar type of osciUations in a strongly magnetised neutron star, 
composed of a fluid core and a solid crust and permeated by a magnetic field with both a toroidal and a poloidal component. 



This mixed magnetic fleld configuration is though to be stable (see Braithwaite & Nordlund (20061, Lasky et al. (20111, 



et aL 



(20111 and references therein). FoUowing the work by (Sotani & Kokkotas 20091 (for the polar part) and 



Ciolfi 



Colaiuda 



fc Kokkotas| ( 2011 1 (for the axial part), we derive a coupled system of equations that involves both the polar perturbation 
functions W and V and the axial function y. The equations are coupled through the toroidal component of the magnetic 
field. After solving those equations with the appropriate boundary conditions and various sets of initial data, we obtain the 
axial and polar spectrum. 

However, we found that both the spectra are strongly influenced by the coupling: the axial spectrum, that in absence of 
couphng had a discrete nature with the lower frequency located around a few hundred Hertz, preserves its discrete nature 
but some very low frequencies appear in the 20-300 Hz range. Those lower frequencies seems to be strictly related to the 
energy transfer between the axial and the polar osciUations via toroidal magnetic field: the lower frequencies found in the 
axial spectrum are the same found in the polar one, with a tiny shift of 1 or 2 Hz. However, the most relevant change is the 
one that concerns the axial spectrum: examining the FFT ampUtude at various points in the star by varying the phase, we 
observe that the axial spectrum has a discrete nature. In fact, we can observe only the edges of the continuum that we found 



Colaiuda & Kokkotas (2011 1: those edges are stiU Alfven modes, in fact they scale with the strength of the magnetic field 



and their restore force is the magnetic force. Due to the presence of the crust tn our model, also 'crustal modes ' appear in 
the spectrum. In the axial case, they are not anymore absorbed by the continuum so they are long-living modes. The only 
case in which a damping of the crustal modes can be observed is if they are in resonance with the edges of the continuum. In 
this case the energy leaks from the crust to the core with a velocity that depends on the strength of the magnetic field. The 
efficiency of the process depends on the magnetic field strength. FinaUy, we made an attempt to resolve the puzzle of the 625 
Hz QPO: we find that in our model, although the spectrum is discrete, at higher frequencies it becomes too dense to allow 
an unambiguous interpretation of this high frequency mode. 

In the future we would like to investigate the effect of superconductivity on the oscillations spectrum as well as non- 
axisymmetric oscillations. Note that a first study of non-axisymmetric oscillations in a magnetised neutron star was made by 
[Lander et al.| ( |2010[ ) , [Lander fc Jones] ( |2011[ ) (with a Newtonian model) and by [Sotani fc Kokkotas| ( |2011[ ) (with a general 
relativistic model): they both found that the spectrum of non-axisymmetric oscillations has also a discrete nature. 



ACKNOWLEDGEMENTS 

We thank Y. Levin, N. Stergioulas and M. Gabler for fruitful discussions and E. Gaertig for providing us his mode recycling 
routine. This work was supported by the German Science Council (DFG) via SFB/TR7. 



APPENDIX: THE SHEAR STRESS TENSOR 

The linearized shear stress tensor SS^i, is defined as: 

= -2/155^" (52) 

and it can be calculated from the formula: 

= e-*55";" (53) 

Here ct^^ is defined as: 

0"mi^ = ^('^M;Q-P° + U^-aP^) - ^Ptii^u'^13 (54) 

with P^i/ = gfj,i, + UfiUv The linearized perturbation of a^^ leads to 

^CTp^ = ^{5u^-aPu + U^.aSP" + SUy-aP^ + Uy:aSP") — ^^Pp^uf^ — ^^P^ivSu^.^ (55) 
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The non-zero components of the tensor 5(7^^ 1 55 1 are: 



5ar0 

Sar^ 
Saee 
5ae4, 
Sasrh 



3r2 



2e-^W.tr - 6 — W,t + Vte + 

r sm H 

2 e"^ 
V,tr --V,t- -^W,te 



-r e amo Y,tr , 



3e- 



-W,t + 2V,te 



-Xt 



1 2 -S . cfl\r 

-r e smB Y,te , 
„-* 

sm ( 

3 



e '^W'.t. + W,t + Vte - 2— 

r sm 



(56) 
(57) 
(58) 
(59) 
(60) 
(61) 



Using the above equations, we can easily calculate the shear stress tensor iJT^" from equation (52 I. Once the shear stress 
tensor is known, we can calculate the shear energy contribution present in equation ( 12 1 for the polar part (i.e. for fj, — r and 
^ = 6). These contributions are given by: 



3 g2#+3A^3 



g2#+A^4 



W.ee + 77 



3 g24>+3A^2 



$' + A' + - VF. - 4 



g2#+2Aj.2 



Ve 



3 g24>+2A^2 



'^'--]V,e + - 
r 



1 /i cot S , , 2/1 cot 6 

y r + 77 



4^ „ , 4^ , ^MCot^ e-2(*+'^)/. 



3g2*j.2 
2g-2(*+A) 



3e2*r4 



3 r4 



3 g2#+2Aj.2 ' 3 g2#+2Aj,2 
$'+A'Vr + 



.^/ A/ 1 2e 
•1' + A' + - - 



3r 



4Me-^^'cotg ^^ 



3r4g2*+A 



W'.se + 



- -]We 



(62) 



(63) 
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